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•  Potential  markets  for  ethanol-fueled  solid  oxide  fuel  cell  power  systems  in  Brazil  are  analyzed. 

•  A  simulation  model  of  this  kind  of  power  system  is  formulated. 

•  The  model  is  run  with  different  system  configurations  to  identify  the  optimal  one. 
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Solid  oxide  fuel  cells  are  globally  recognized  as  a  very  promising  technology  in  the  area  of  highly  efficient 
electricity  generation  with  a  low  environmental  impact.  This  technology  can  be  advantageously  imple¬ 
mented  in  many  situations  in  Brazil  and  it  is  well  suited  to  the  use  of  ethanol  as  a  primary  energy  source, 
an  important  feature  given  the  highly  developed  Brazilian  ethanol  industry.  In  this  perspective,  a 
simplified  mathematical  model  is  developed  for  a  fuel  cell  and  its  balance  of  plant,  in  order  to  identify  the 
optimal  system  structure  and  the  most  convenient  values  for  the  operational  parameters,  with  the  aim  of 
maximizing  the  global  electric  efficiency.  In  this  way  it  is  discovered  the  best  operational  configuration 
for  the  desired  application,  which  is  the  distributed  generation  in  the  concession  area  of  the  electricity 
distribution  company  Elektro.  The  data  regarding  this  configuration  are  required  for  the  continuation  of 
the  research  project,  i.e.  the  development  of  a  prototype,  a  cost  analysis  of  the  developed  system  and  a 
detailed  perspective  of  the  market  opportunities  in  Brazil. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  field  of  new  energy  converters  has  recently  been  the  target 
of  intense  research  activities,  with  the  main  goal  of  improving 
systems'  energy  efficiency,  while  facing  the  challenges  related 
pollutant  emissions  reduction.  In  this  scenario,  one  of  the  most 
promising  technologies  is  the  Fuel  Cell  (FC).  Fuel  cells  are  electric 
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energy  generators  characterized  by  their  high  efficiency  (about  50% 
based  on  the  net  energy  produced  [1])  and  very  low  pollutant 
emissions.  The  electric  energy  generation  results  from  direct  con¬ 
version  of  the  fuel's  chemical  energy  through  electrochemical  re¬ 
actions.  In  addition  to  the  fuel  cells  themselves,  arranged  in  stacks, 
electric  generation  systems  based  on  FCs  require  various  auxiliary 
devices,  collectively  called  Balance  of  Plant  (BoP)  of  the  fuel  cell 
stack. 

One  type  of  FC  that  attracts  great  interest  in  the  scientific 
community  is  the  Solid  Oxide  Fuel  Cell  (SOFC),  which  uses  a  solid 
ceramic  electrolyte  capable  of  conducting  oxygen  ions  (02~)  from 
the  cathode  to  the  anode  [1],  The  operating  temperature  varies 
between  650  °C  and  1000  °C.  The  fuel  can  be  hydrogen,  biogas  or 
simple  hydrocarbons,  which  can  decompose  to  CO  and  H2  inside 
the  SOFC  anode,  or  synthesis  gas  obtained  through  previous 
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chemical  treatment  of  complex  hydrocarbons.  Two  important  ad¬ 
vantages  of  this  kind  of  FC  are  the  capacity  to  work  with  fuels  other 
than  pure  hydrogen  and  the  high  temperature  of  the  heat  gener¬ 
ated  as  a  byproduct  of  the  electric  energy  production.  Due  to  their 
peculiar  characteristics,  SOFCs  are  particularly  well  suited  to  use  in 
high  efficiency  stationary  Combined  Heat  and  Power  (CHP)  sys¬ 
tems,  with  power  ranging  from  500  W  up  to  several  MW.  In  fact,  the 
energetic  performance  of  fuel  cell  based  generation  systems  is 
approximately  a  constant  function  of  the  produced  power. 

One  of  the  most  interesting  applications  of  low  power  SOFCs  is 
in  micro-CHPs  (combined  heat  and  electric  power  generators 
characterized  by  nominal  electric  power  lower  than  50  kW,  ac¬ 
cording  to  the  European  Union's  definition  [2]).  This  is  a  techno¬ 
logical  area  where  relevant  development  efforts  are  concentrated, 
particularly  aiming  to  obtain  more  efficient,  sustainable  and  stable 
national  electric  generation  grids,  by  means  of  the  widespread  use 
of  high  efficiency  micro-CHP  units  in  a  distributed  generation 
scheme,  as  evidenced  in  the  literature  [3  .  However,  to  reach  a  true 
beneficial  effect,  it  is  vitally  important  that  these  micro-CHP  units 
have  an  efficiency  similar  to  that  observed  in  modern  centralized 
generation  systems.  This  requirement  cannot  be  easily  satisfied 
using  conventional  technologies,  but  it  can  be  reached  using  FC 
systems,  which  present  high  efficiencies  even  when  working  at  low 
power  levels.  As  a  consequence,  fuel  cell  systems  are  potential 
candidates  for  massive  diffusion  of  micro-CHP  systems.  In  this  area, 
SOFC  technology  seems  to  be  particularly  promising  [4],  offering 
the  previously  mentioned  advantages  over  other  types  of  FCs. 
During  the  last  ten  years,  the  research  in  this  field  has  been  intense, 
leading  to  the  creation  of  several  prototypes  and  pre-commercial 
products.  Several  installation  programs  were  initiated  in  order  to 
demonstrate  the  technology  [5,6],  and  some  companies  are  in  the 
first  commercialization  stage  of  SOFC-based  methane-fueled 
micro-CHPs,  introduced  in  the  market  as  substitutes  for  the  tradi¬ 
tional  domestic  heaters,  widely  used  for  residential  applications. 
Given  the  proposed  application,  these  devices  are  usually  designed 
in  a  way  to  maximize  global  efficiency  (sum  of  thermal  and  elec¬ 
trical  efficiencies)  and  to  assure  enough  thermal  power,  treating  the 
electric  energy  just  like  as  a  useful  byproduct  of  heat  production. 

2.  Potential  markets  in  Brazil 

In  Brazil,  energy  generation  and  its  use  have  peculiar  profiles, 
very  different  from  the  ones  that  characterize  most  other 
economically  important  countries  on  the  global  scale.  One  impor¬ 
tant  characteristic  is  the  high  percentage  of  electric  energy  that  is 
produced  through  renewable  sources  (89%),  mainly  hydropower 
and  biomass,  with  a  small  wind  power  contribution  [7  .  Another 
important  aspect  is  the  minimal  use  of  energy  for  residential 
heating,  limited  to  small  areas  of  the  country  and  only  for  few  days 
a  year.  Moreover,  the  country  presents  the  need  for  a  considerable 
expansion  of  the  energy  production,  transmission  and  distribution 
infrastructure,  to  support  the  country's  economic  growth  and 
extend  the  availability  of  electrical  power  in  remote  areas.  To  tackle 
this  last  issue,  a  federal  program  was  launched:  between  2003  and 
2013  more  than  3  million  new  electric  connections  were  imple¬ 
mented,  with  an  investment  of  R$  20  billion  [8],  Despite  this  in¬ 
vestment,  some  areas  of  the  country  still  lack  a  connection  to  the 
national  transmission  grid,  and  a  substantial  growth  of  demand  is 
expected.  Therefore,  it  is  probable  that  the  country's  power  supply 
needs  will  grow  even  more  in  the  next  years,  requiring  additional 
grid  reinforcement  or  alternative  solutions  such  as  distributed 
generation.  Finally,  Brazil  possesses  an  important  bioethanol  pro¬ 
duction  capacity  and  distribution  infrastructure  [7],  well  estab¬ 
lished  in  every  Brazilian  state,  including  both  large  industries  and 
small  family-operated  agro-industries.  The  primary  energy 


generation  from  sugar  cane  products  showed  steep  growth  in  the 
last  years  and  in  2009  surpassed  the  hydroelectric  energy  genera¬ 
tion  in  the  country. 

Given  this  scenario,  the  use  of  SOFC  micro-CHP  systems  would 
also  present  some  peculiarities,  such  as  the  use  of  ethanol  as  fuel, 
instead  of  methane,  and  the  operation  essentially  dedicated  to 
electric  production,  even  though  some  areas  of  the  country  present 
periods  of  cold  weather,  when  the  residential  cogeneration  of 
electricity  and  heat  could  be  interesting.  There  are  many  situations 
in  Brazil  where  the  use  of  SOFC  systems  could  bring  environmental 
and  economic  benefits.  A  first  application  is  the  substitution  of 
generators  characterized  by  high  environmental  impact  (in  most 
cases  internal  combustion  engines).  These  are  currently  used  in 
some  specific  cases,  including  the  generation  of  power  when  the 
electric  energy  reaches  its  peak  in  demand  (and  as  a  consequence 
the  cost  is  higher),  or  the  supply  of  electricity  to  remote  places  not 
connected  to  the  national  transmission  grid,  or  even  the  strength¬ 
ening  of  load  covering  capacities  in  small  groups  of  grid  users 
distant  from  the  main  company  grid  that  do  not  present  a  regular 
consumption  high  enough  to  justify  the  installation  of  new  higher 
capacity  cables.  It  would  be  possible  to  replace  the  above 
mentioned  generators  with  SOFC  systems,  maintaining  the  same 
service  and  infrastructure  quality,  but  markedly  increasing  gener¬ 
ation  efficiency  and  at  the  same  time  reducing  pollutant  emissions. 
The  Brazilian  market  for  peak-hour  additional  generation  capacity 
is  estimated  to  lie  between  3  and  10  GW  [9],  Until  recently,  the 
energy  generated  to  meet  this  peak-hour  demand  was  produced 
through  the  use  of  internal  combustion  engines,  but  an  increase  in 
the  price  of  diesel  and  a  drop  in  the  electricity  peak  price  have 
made  this  solution  less  economically  viable.  On  the  other  hand, 
ethanol-fueled  SOFC  systems  could  offer  better  economic  perfor¬ 
mance,  due  to  this  type  of  equipment's  higher  electric  efficiency. 
Furthermore,  the  electric  energy  demand  projected  for  the  next 
years  can  be  used  to  evaluate  the  potential  market  size  for 
distributed  generation.  The  expected  increase  in  Brazilian  elec¬ 
tricity  demand  between  2012  and  2022  has  been  estimated  around 
275  TWh,  including  grid  supplied  power  and  autoproduction  (in¬ 
dustrial  cogeneration  and  distributed  generation)  [10].  This  value 
represents  an  increase  of  approximately  55%  whit  respect  to  the 
current  figures.  Admitting  that  a  large  part  of  this  considerable 
demand  increase  will  be  satisfied  using  conventional  methods 
(namely,  expansion  of  centralized  generation  and  power  trans¬ 
mission  grids),  the  potential  market  for  distributed  generation  in 
Brazil  is  still  very  important,  especially  considering  that  the  size  of 
the  country  implies  high  costs  of  transmission  lines.  For  example, 
the  current  cost  of  a  345  kV  line  can  be  evaluated  around  584,000 
R$  km-1  [11  .Asa  consequence,  in  several  cases  the  cost  of  the  lines 
required  to  connect  remote  locations  to  the  national  transmission 
grid,  or  to  increase  the  power  capacity  of  existing  connections,  can 
make  the  dissemination  of  distributed  generators  with  higher  cost 
per  installed  kW  economically  viable.  In  this  scenario,  ethanol- 
fueled  SOFC  micro-CHPs  could  offer  unique  advantages  over  other 
technologies,  namely  the  possibility  to  generate  energy  on  demand 
using  a  renewable  energy  source.  In  addition  to  what  was  stated  up 
to  now  in  this  section,  it  should  be  noted  that,  even  though 
centralized  electric  generation  in  Brazil  basically  relies  on  renew¬ 
able  sources,  its  environmental  impact  is  not  negligible.  Being 
mostly  of  hydro-electrical  origin,  it  involves  a  considerable  pro¬ 
duction  of  methane  and  carbon  dioxide  from  organic  material 
decomposition  in  water  reservoirs  [12];  in  fact,  the  extension  of 
water  accumulation  lakes  can  be  used  to  quantify  the  environ¬ 
mental  impact  of  these  plants  [13],  Moreover,  since  it  is  centralized, 
this  kind  of  generation  requires  a  big  transmission  and  distribution 
infrastructure  that  presents  high  financial  and  environmental  costs, 
and  creates  relevant  losses  in  efficiency.  As  a  consequence,  the  use 
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of  SOFC-based  micro-CHPs  fueled  by  bioethanol  in  a  distributed 
generation  pattern  supplying  base  power  to  the  grid  could  repre¬ 
sent  a  practical  way  of  generating  additional  renewable  electricity 
with  a  reduced  environmental  impact. 

Given  the  situation  described,  a  research  project  was  initiated, 
aiming  at  the  development  of  a  hybrid  system  based  on  an  ethanol- 
fueled  SOFC  and  with  the  main  objective  of  generating  electric 
energy  in  the  micro-CHP  power  range.  The  term  “hybrid”  means 
that  the  SOFC  should  produce  a  constant  net  power  of  500  W,  as 
part  of  a  more  complex  system  that  also  includes  a  battery  pack 
used  to  compensate  for  the  electric  load  variations,  thereby  main¬ 
taining  the  fuel  cell  power  output  essentially  constant.  In  this  way, 
it  is  possible  to  reduce  the  power  transients  the  fuel  cell  needs  to 
sustain,  increasing  its  life  span  14],  In  fact,  frequent  variations  in 
the  fuel  cells’  operational  conditions  can  lead  to  accelerate  perfor¬ 
mance  degradation  [15  .  Another  important  characteristic  of  the 
system  will  be  the  direct  use  of  ethanol,  without  previous  chemical 
processing.  This  will  be  possible  with  the  use  of  a  recently  devel¬ 
oped  anode  material  for  SOFC  fabrication  [16,17  which  possesses 
electrocatalytic  action  for  the  direct  oxidation  of  ethanol  without 
solid  carbon  formation  inside  the  fuel  cells.  Also  the  BoP  structure, 
materials  and  management  strategy  contribute  in  avoiding  carbon 
deposition  problems.  This  characteristic  will  be  an  important  dif¬ 
ference  with  respect  to  other  existing  prototypes,  where  the  fuel 
goes  through  a  series  of  chemical  treatments  before  entering  the 
fuel  cell  stack.  The  application  of  these  treatments  would  imply  the 
use  of  one  or  more  additional  reactors  that  increase  the  cost  and  the 
complexity  of  the  system. 

The  first  step  of  the  research  consisted  in  conceiving  a  simplified 
mathematical  model  of  the  SOFC-BoP  ensemble  (excluding  the 
batteries  and  the  interface  generator/grid)  and  using  it  to  evaluate 
some  possible  system  construction  topologies,  looking  for  the 
optimal  values  of  the  operational  parameters  in  each  different  case. 
The  data  obtained  will  be  used  in  the  further  steps  of  the  research 
to  fix  the  specifications  of  the  system  and  to  build  an  experimental 
prototype.  The  applied  model  is  purely  thermodynamic,  not  con¬ 
taining  any  economic  consideration  about  the  cost  of  the  system. 

3.  Mathematical  model 

The  system's  mathematical  model  includes  the  main  compo¬ 
nents  that  have  been  proposed  for  use  in  the  construction  of  the 
SOFC-BoP  ensemble  prototype.  The  list  of  components  was  created 
through  basic  project  activity  of  the  system  and  comparison  with 
some  structures  described  in  the  literature  ([18,19]).  Some  of  the 
chosen  components  (labeled  as  “optional”)  are  not  present  in  all  the 
versions  that  were  considered  as  alternatives  for  system 
fabrication: 

•  ethanol  pump,  used  to  feed  fuel  to  the  stack; 

•  ethanol  tank; 

•  water  pump  (optional),  can  be  employed  to  create  a  water  flow, 
used  to  prevent  the  formation  and  the  deposition  of  solid  carbon 
in  the  system's  pipes  and  inside  the  SOFC; 

•  water  tank  (optional); 

•  blower,  used  to  generate  the  air  flow  needed  for  fuel  cell 
reactions; 

•  after-burner,  a  combustor  where  the  excess  fuel  which  did  not 
react  in  the  SOFC  burns  completely; 

•  fuel  and  air  pre-heating  heat  exchangers; 

•  recirculation  ejectors  (optional),  used  to  promote  a  partial 
recirculation  of  anode  and  cathode  outlet  gases  to  the  anodic 
and  cathodic  inlets  respectively.  In  the  ejectors,  the  fuel  (or  air) 
flow  created  by  the  ethanol  pump  (or  the  air  blower)  transfers 
part  of  its  momentum  to  a  portion  of  the  flow  that  exits  from  the 


anode  (or  the  cathode),  generating  the  recirculation  of  a  portion 
of  the  gas,  which  reenters  the  SOFC; 

•  power  electronic  device,  properly  conditions  the  electric  power 
produced  by  the  fuel  cell; 

•  heat  recovery  heat  exchanger,  uses  the  residual  energy  content 
of  after-burner  outlet  gases  to  heat  a  flow  of  water. 

Fig.  1  shows  the  basic  layout;  numeric  indeces  are  used  in  the 
model  description  to  indicate  different  system  points. 

The  modeled  generation  system  has  the  main  goal  of  producing 
electric  energy  at  constant  power,  while  a  battery  bank,  not 
included  in  the  model,  compensates  for  load  power  fluctuations.  As 
a  consequence,  in  this  stage  of  the  development  the  model  de¬ 
scribes  only  the  stationary  system  regime. 

3.1.  SOFC  model 

The  electricity  generation  device  of  the  system  is  a  SOFC  stack, 
an  electrochemical  equipment  capable  of  converting  ethanol's 
chemical  energy  directly  into  electric  energy.  Actually,  ethanol  is 
decomposed  in  simpler  compounds  through  different  reactions 
before  it  reaches  the  FC  anode.  The  result  of  these  reactions  is  the 
production  of  hydrogen,  water,  carbon  monoxide  and  carbon  di¬ 
oxide.  The  elements  that  take  part  in  the  electrochemical  reactions 
are  essentially  carbon  monoxide  and  hydrogen  in  the  anode,  (1) 
and  (2),  and  oxygen  in  the  cathode,  (3): 


h2+o2-^ 

►  H20  +  2e~ 

(1) 

CO  +  O2- 

-+C02  +  2e~ 

(2) 

^02  +  2e- 

-+02- 

(3) 

In  this  model,  it  is  assumed  that  all  ethanol  is  converted  to  CO 
and  H2  inside  fuel  cells  anodic  compartments  or  in  gas  feeding 
pipes,  before  reaching  the  active  sites  where  the  chemical  reactions 
take  place.  This  result  can  be  obtained  in  a  real  system  through 
different  strategies  [20]: 

•  external  reforming  —  the  fuel,  together  with  steam,  goes 
through  a  reactor  that  contains  a  catalyst  capable  of  promoting 
reactions  that  break  down  the  ethanol  molecule.  Common  re¬ 
actions  are  called  steam  reforming  (4)  and  water-gas  shift  (5): 


C2  H60  +  H2  0  ^  2CO  +  4H2 

(4) 

co  +  h2o^co2  +  h2 

(5) 

The  first  reaction  requires  high  temperatures  and  is  strongly 
endothermic;  the  second  is  mildly  exothermic.  This  means  that 
generally  it  is  necessary  to  supply  heat  to  the  reactor  in  order  to 
maintain  its  thermal  equilibrium.  Hence,  the  presence  of  this 
external  reactor  considerably  complicates  the  system’s  construc¬ 
tion  and  control,  increasing  its  cost.  Due  to  these  considerations, 
the  described  project  focuses  on  the  research  of  solutions  that  do 
not  involve  external  reforming.  Therefore,  this  option  was  not 
included  in  the  alternatives  explored  by  the  model; 

•  internal  reforming  —  ethanol  is  fed  directly  to  the  SOFC  blended 
with  steam.  In  this  case  there  is  no  dedicated  reforming  reactor: 
(4)  and  (5)  happen  in  the  anodic  compartments.  In  this  way,  it  is 
possible  to  obtain  an  optimal  thermal  balance  between  the  heat 
produced  by  the  electrochemical  reactions  (1)— (3),  exothermic 


A.  Coralli  et  al.  /  Journal  of  Power  Sources  269  (2014)  632-644 


635 


Fig.  1.  Basic  layout  of  SOFC-BoP  ensemble.  In  the  scheme,  all  the  components  listed  in  this  section  are  represented,  together  with  the  connections  between  them.  Dashed  lines 
indicate  optional  piping  connections. 


reactions,  and  the  heat  required  by  the  steam  reforming  (4).  The 
resulting  system  is  simpler  and  cheaper  than  in  the  previous 
case,  but  it  still  needs  a  steam  source.  The  viability  of  this  option 
was  explored  in  the  literature  [21]; 

•  direct  utilization  —  the  SOFC  stack  is  fed  with  pure  ethanol, 
without  any  steam  addition.  However,  the  ethanol  molecule  is 
ripped  apart  by  pyrolysis  reactions  before  reaching  the  interior  of 
the  porous  matrix  of  the  anode,  producing  simpler  compounds 
that  can  react  at  the  active  sites  that  exist  in  the  region  of  the  anode 
where  the  electrochemical  reaction  takes  place.  When  the  stack 
reaches  the  steady  state,  water  produced  by  the  oxidation  of 
hydrogen  (1)  is  present  at  the  anode.  Therefore,  in  this  condition, 
reactions  (4)  and  (5)  can  happen,  even  though  no  steam  is  being 
added  from  outside  the  stack  (as  verified  in  Ref.  [22]).  This  solution 
guarantees  lower  system  complexity.  On  the  other  hand,  when  a 
SOFC  stack  is  operated  in  this  way  there  is  generally  a  possibility  of 
solid  carbon  formation  inside  the  stack  and  the  high  temperature 
pipes.  The  components  of  the  system  which  are  potentially 
affected  by  soot  formation  are  the  parts  of  the  fuel  inlet  pipes 
where  the  temperature  is  high  enough  to  favor  fuel  pyrolysis  re¬ 
actions,  the  fuel  inlet  manifold  of  the  stack  and  the  anode  com¬ 
partments  of  the  fuel  cells.  Moreover,  pyrolysis  reactions  can  be 
favored  by  the  contact  between  the  fuel  and  some  elements  or 
compounds,  which  can  act  as  catalysts  and  lower  the  temperature 
range  in  which  solid  carbon  formation  is  considerable.  As  a 
consequence,  a  possible  way  to  eliminate  this  problem  is  the  use  of 
materials  that  do  not  favor  or  avoid  carbon  deposition.  In  the 
present  paper,  it  is  assumed  that  the  use  of  an  innovative  material 
([16,  [17]),  already  mentioned  in  Section  2,  and  the  BoP  structure, 
materials  and  management  strategy  make  it  possible  to  obtain  the 
optimal  situation  in  which  there  is  no  carbon  deposition,  neither 
in  the  stack  nor  in  other  parts  of  the  system. 

The  present  model  considers  only  the  system's  steady  state 
condition.  In  this  situation,  the  SOFC  anode  contains  water  also  in 
the  case  of  the  direct  utilization  of  ethanol.  For  simplicity,  the 
model  is  applied  only  with  the  combinations  of  operational  pa¬ 
rameters  that  guarantee  that  the  amount  of  water  in  the  stack 
anodes  is  enough  to  convert  through  reaction  (4)  all  of  the  ethanol 


supplied.  In  addition  to  that,  in  normal  operating  conditions  reac¬ 
tion  (5)  happens  approximately  four  times  faster  than  (2)  [23).  As  a 
consequence,  in  this  model  the  water-gas  shift  reaction  is 
considered  as  the  only  mechanism  capable  of  consuming  carbon 
monoxide  (electrochemical  oxidation  is  neglected).  The  global 
result  of  (1)— (5)  is  therefore  ethanol  conversion  to  CO2  and  H2  with 
subsequent  electrochemical  oxidation  of  the  hydrogen  produced. 
The  magnitude  of  the  errors  introduced  with  these  hypotheses  will 
be  quantified  in  the  subsequent  steps  of  the  research. 

The  model  is  built  assuming  that  (1)— (5)  happen  at  the  same 
time  and  that  all  carbon  contained  in  the  ethanol  reacts,  forming 
carbon  dioxide.  In  this  way,  it  is  possible  to  write  global  reactions, 
also  considering  the  possibility  of  partial  recirculation  of  reaction 
products  between  anode  and  cathode  inlets  and  outlets.  The  for¬ 
mulas  were  obtained  applying  (4)  and  (5)  for  all  ethanol  supplied, 
and  (1)— (3)  only  for  the  amount  allowed  by  transferred  oxygen 
ions.  Both  anodic  (6)  and  cathodic  (7)  reactions  were  written  using 
as  calculation  basis  six  of  the  oxygen  ions  among  those  transferred 
from  the  cathode  to  the  anode  through  the  electrolyte: 

a  C2H6O  ■  ^  H2O  +  /?C02  +  5H2  +  6O2 

-  3a  +  6^) H20  +  {P  +  2)C02  +  (6a  +  5  -  6)H2  +  12e 

(6) 

/j,C>2  +  vN2  +  12e~  — *  (/j.  —  3)02  +  FN2  +  602~  (7) 

where  S/£  is  the  ratio  between  the  number  of  moles  of  steam  and 
ethanol  entering  the  stack,  a,  f),  5,  ji  and  v  are  the  stoichiometric 
coefficients  of  ethanol,  carbon  dioxide,  hydrogen,  oxygen  and  ni¬ 
trogen,  respectively.  To  calculate  these  coefficients,  a  balance  of 
available  electrons  at  the  anodic  inlet  was  used: 

12a  +  2<5  =  127fuel  (8) 

7 fuei  indicates  the  ratio  between  the  amount  of  fuel  supplied  and 
the  stoichiometric  amount,  which  is  the  minimal  amount  of  fuel 
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required  to  complete  the  reactions  with  the  oxygen  ions.  Afuei  is  the 
inverse  of  the  fuel  utilization  rate,  another  parameter  widely  used 
in  the  literature.  An  additional  relation  is  introduced  to  allow  for 
the  calculation  of  the  stoichiometric  coefficients  and  is  written 
considering  that  the  anodic  recycling  fraction  /rec.a.  defined  as  the 
ratio  between  the  amount  of  a  chemical  compound  entering  the 
anode  as  a  result  of  recirculation  and  the  amount  of  the  same 
compound  leaving  the  anode,  is  the  same  for  all  compounds  that 
are  not  introduced  from  outside  the  stack: 


Vc  =  EN  -  i-Rc  -  a  -  b  ln(i)  (14) 

with  Vc  single  cell  voltage,  i  current  density  in  A  citT2,  Rc  ohmic 
resistance,  a  and  b  constants  related  to  the  activation  losses.  Rc,  a 
and  b  are  calculated  as  a  function  of  temperature  through  empirical 
expressions  found  in  the  literature  for  the  case  of  a  typical  solid 
oxide  cell  [23],  If  Nc  and  Ac  are  respectively  the  number  of  cells  in 
the  stack  and  the  active  area  of  a  single  cell,  the  electric  power 
produced  can  be  expressed  by: 


/reca  = 


(3  _  -/H20,e) 


6a +  5  —  6  13  +  2 


4  -3a  +  6 


(9) 


This  is  not  valid  for  the  substances  that  are  introduced  from 
outside  the  stack  (ethanol,  water,  oxygen  and  nitrogen),  since  in 
this  case  it  is  necessary  to  include  also  the  external  flow  in  the  mass 
balance.  Two  additional  project  parameters  are  used  to  determine 
the  values  related  to  these  compounds  as  well:  the  fraction  of  water 
flow  introduced  from  outside  the  stack,  /h2o.  and  the  cathodic 
recirculation  fraction /rec.c-  This  last  parameter,  defined  in  the  same 
way  as /rec,a.  can  also  be  used  to  evaluate  the  fraction  of  oxygen  that 
is  introduced  from  outside  the  stack: 


fo2,e 


rcp2,e 

^02,in 


f REC,C 


(10) 


which  is  the  ratio  between  the  molar  flow  rate  of  oxygen  intro¬ 
duced  in  the  stack  from  outside  n0i  e  and  the  total  molar  flow  rate 
of  oxygen  entering  the  cathode  compartments  n0l  in.  In  addition,  it 
is  assumed  that  the  cathodic  fluid  is  air  (this  implies  that  approx¬ 
imately  3.76  mol  of  nitrogen  are  introduced  in  the  stack  for  every 
mole  of  oxygen).  Finally,  it  is  possible  to  introduce  the  parameter 
Aq2,  analogous  to  Afuei,  which  indicates  the  ratio  between  the 
amount  of  oxygen  actually  supplied  and  the  stoichiometric  amount, 
which  is  the  minimal  amount  of  oxygen  required  to  complete  the 
reactions  with  the  electrons  coming  from  the  anode  through  the 
electrolyte.  From  reaction  (7),  it  can  be  concluded  that  Aq2  can  be 
expressed  as  in  equation  (11): 


^=1  (11) 

The  stoichiometric  coefficients  introduced  in  (6)  and  (7)  give  the 
molar  fractions  of  all  chemical  substances  contained  in  the  gas  flow 
exiting  from  the  stack,  according  to  following  relation: 


l*ei,S0FC  —  VI  =  NcVciAc  (15) 

Vc  and  i  are  calculated  fixing  Nc,  Ac  and  Pei  as  project  data.  The 
required  molar  flow  rate  of  ethanol  is  related  to  the  total  current  / 
(product  of  current  density  and  cell's  active  area)  through  the 
following  equation  that  is  obtained  considering  the  electronic 
balance  (8)  and  the  reaction  (6): 


I  Nc- a 
nc2H6o  =  -J2J+- 


(16) 


The  molar  flow  rates  of  the  other  compounds  present  in  the 
anode  inlet  and  outlet  flows  are  calculated  according  to  the  stoi¬ 
chiometry  indicated  in  (6),  using  the  generic  relation  (17): 


Hi  =  ric2H6cr 


Q 

cc2h6o 


ri  Ci 
nC2H60'- 


(17) 


where  Cc2h6o  is  the  stoichiometric  coefficient  of  ethanol  (equal  to  a) 
and  Ci  is  the  stoichiometric  coefficient  of  the  chemical  specie  i, 
related  to  the  inlet  or  outlet  molar  flow  rate  n\  of  the  chemical 
specie  i.  In  a  similar  way,  the  molar  flow  rate  of  oxygen  introduced 
in  the  stack  from  outside  can  be  calculated  using  equation  (18), 
analogous  to  equation  (16): 


ho2,e  =  [a°2  -/r£C,C'  (Ao2  -  1)]  (18) 

and  the  molar  flow  rate  of  nitrogen  introduced  in  the  stack  from 
outside  is  calculated  in  equation  (19)  taking  into  account  the  defi¬ 
nition  of fo2  e  expressed  by  equation  (10): 


»iN2,e  =3.76  ri02e  =  3.76  /o2,e-ho2  in  (19) 

The  total  chemical  power  introduced  in  the  SOFC  through  the 
fuel  flow  results  from  the  relation: 


*f=A  (12) 

2-4  J 

C,  and  Cj  are  the  stoichiometric  coefficients  of  compounds  i  and  j, 
and  the  sum  is  extended  to  all  compounds  that  are  present  in  the 
outlet  flow.  After  determining  the  reactions’  stoichiometry,  the 
Nernst  equation  was  used  in  order  to  define  the  theoretical  open 
circuit  voltage  for  each  single  cell: 


£n 


AG° 

n-F 


RTS 

n-F 


■In 


(13) 


An  explanation  of  each  term  can  be  found  in  common  fuel  cell 
literature  [1  ].  When  a  load  is  applied,  this  maximum  voltage  level  is 
reduced  due  to  ohmic  and  activation  losses,  as  indicated  in  equa¬ 
tion  (14)  (concentration  losses  are  neglected,  assuming  that  the 
stack  does  not  work  close  to  maximum  current  condition): 


ftot  =  ~AHC2h6o  riC2Hf.0  (20) 

where  AHc2h6o  is  ethanol’s  enthalpy  of  combustion.  The  stack's 
electric  efficiency  is: 


^el.SOFC 

?,eLS0FC=  Ptot 


(21) 


As  previously  stated,  in  the  stack  each  mole  of  ethanol  reacts 
with  steam  according  to  (4)  and  (5),  and  part  of  the  produced 
hydrogen  reacts  in  (1).  It  is  also  necessary  to  include  in  the  model 
the  effects  of  the  reactions  over  SOFC's  thermal  balance.  By  intro¬ 
ducing  hydrogen's  enthalpy  of  combustion  A HH2,  the  sum  of  the 
enthalpies  of  reaction  associated  with  steam  reforming  and  water 
gas  shift  reactions  A Hsr,  and  the  thermal  power  absorbed  by  the 
gases  circulating  in  the  cells  through  convection  Qconv.  it  is  possible 
to  obtain  the  thermal  power  that  is  produced  by  the  stack  using  the 
following  energy  balance: 
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QsOFC  =  -  AWc2H6o  ric2H60  +  AWh2  '  ^H2  -  Pel  -  AHsr '  'iQHeO  =  Qconv 

(22) 

where  nHi  is  hydrogen  flow  at  SOFC  outlet.  If  the  thermal  powers 
associated  with  anode  and  cathode  inlet  gas  flows  are  indicated 
with  Q2  and  Q4,  the  sum  of  the  outlet  flows  thermal  powers  is: 

05  =  02  +  0.4  +  Qconv  (23) 

An  estimate  of  the  convective  thermal  transfer  coefficients  in 
SOFC  cells  channels  [24]  shows  that  the  thermal  exchange  is  effi¬ 
cient  enough  to  increase  the  temperature  of  the  outlet  gases  to  the 
same  level  as  that  of  the  solid  components  of  the  stack.  Inlet 
temperatures  are  project  parameters.  Knowing  all  molar  flows,  it  is 
possible  to  calculate  the  heat  flows  using  the  equation: 

Qk  =  Tk-(Y)cpj(Tk)-nlk')  (24) 

i 

where  k  indicates  the  considered  position,  and  the  sum  is  extended 
to  all  present  chemical  compounds.  cPj(Ti<)  indicates  the  specific 
heat  at  constant  pressure  of  compound  j  at  temperature  T\(. 

3.2.  After-burner  model 

In  order  to  increase  the  thermal  energy  content  of  SOFC  outlet 
gases  and  to  avoid  the  emission  in  the  environment  of  harmful 
substances  like  carbon  monoxide,  the  use  of  an  after-burner  was 
considered.  In  this  component,  placed  after  the  stack,  residual 
inflammable  gases  burn  according  to  the  equation: 

?H20  +JCO2  +  ;H2  +  r02  +  tN2  ->;'N2  +  (£  +  ?)H20  +JC  02 


(25) 

where  £,  q,  r,  e  <p  are  the  stoichiometric  coefficients  of  the 
different  substances.  Formula  (25)  was  written  under  the 
assumption  that  the  oxygen  that  feeds  the  SOFC  is  enough  to  sus¬ 
tain  the  electrochemical  reaction  and  completely  burns  the  portion 
of  the  fuel  that  does  not  react  in  the  cells.  Using  the  stoichiometric 
coefficients  and  knowing  the  composition  of  the  after-burner  inlet 
flow,  it  is  possible  to  calculate  the  molar  flows  of  the  chemical 
compounds  leaving  the  after-burner.  The  total  molar  flow  of  the 
outlet  gases  ns  is  the  sum  of  these  flows.  Since  power  Qsofc  may  be 
a  negative  value  (which  means  that  the  SOFC  requires  thermal 
energy  in  order  to  maintain  its  thermal  equilibrium),  it  is  assumed 
that  the  after-burner  is  in  direct  thermal  contact  with  the  stack, 
exchanging  heat  flow  Qsofc-  In  this  way,  the  SOFC  internal  tem¬ 
perature  is  considered  the  same  for  every  simulated  condition. 
Indicating  with  Qs  and  Qd,AB  the  thermal  power  associated  with  the 
after-burner  outlet  gases  and  the  thermal  power  that  is  dispersed 
through  the  perimeter  of  the  after-burner  respectively,  it  is  possible 
to  formulate  the  following  energy  balance: 

06  =  Os  -  AHh2  •rcH2  -  Od.AB  +  Qsofc  (26) 


3.3.  Blower  and  pumps  models 

The  power  consumed  by  the  fuel  pump  (and  eventually  by  the 
water  pump)  can  be  expressed  as  a  function  of  the  molar  flow  of  the 
fluid,  the  specific  enthalpy  variation  between  pump  inlet  and  outlet 
A hf,  and  the  isentropic  efficiency  of  the  machine  tjp.  Alternatively, 
the  same  quantity  can  be  represented  using  the  volumetric  flow  Vf 
and  the  pressure  variation  A pf  (the  temperature  variation  in  a 
pump  is  negligible): 

Vp  Vp 

An  analogous  equation  can  be  developed  for  the  power 
consumed  by  the  air  blower  [14],  It  is  possible  to  express  the  blower 
outlet  flow  temperature  for  an  isentropic  compression  process 
using  the  following  equation: 


^3, isentropic  —  ^amb  '  (p  (29) 

As  a  consequence,  the  power  consumed  by  the  blower  can  be 
expressed  as: 


Ph  = 


^air'^^b  ^air 'cair(^amb) '^amb 


Vb 


Vb 


P^)"  -1 

,Pamb/ 


(30) 


where  7jb  is  the  blower’s  isentropic  efficiency,  p0  is  the  blower's 
outlet  pressure  and  y  is  the  adiabatic  expansion  coefficient. 

To  execute  the  calculations,  the  isentropic  efficiency  r/b  was 
considered  equal  to  0.7  and  the  pressure  drops  in  the  various 
components  were  fixed  according  to  Table  1. 

3.4.  Heat  exchangers  model 

The  basic  layout  of  the  system  includes  three  heat  exchangers, 
used  to  recover  part  of  the  energy  content  of  the  outlet  gases,  pre¬ 
heating  the  fluids  that  are  fed  to  the  stack  and  heating  a  water  flow, 
which  can  be  used  by  an  external  user.  The  following  energy  bal¬ 
ance  can  be  applied  to  each  heat  exchanger: 

Qh.i  ~  Qh.o  +  Qh.vap  =  Qc.o  —  Qc.i  +  0d,TR  +  Qcvap  (31 ) 

In  equation  (31 ),  the  heat  flows  related  to  the  hot  fluid  and  to  the 
cold  fluid  are  indicated,  respectively,  with  subscripts  h  and  c,  while 
subscripts  i  and  o  indicate  the  heat  exchanger's  inlet  and  outlet 
flows.  Qd.TR  is  the  thermal  power  dispersed  through  the  heat  ex¬ 
changer's  insulation.  Qh.vap  and  Qc.vap  are  the  latent  heats  of 
vaporization  of  hot  and  cold  fluid,  thermal  powers  that  have  to  be 
removed  or  supplied  in  order  to  condensate  the  hot  fluid  or 
vaporize  the  cold  fluid.  Applying  this  balance  with  equation  (24),  it 
is  possible  to  calculate  inlet  and  outlet  temperatures  for  all  three 
heat  exchangers,  as  well  as  the  thermal  power  absorbed  by  the 
water  in  the  heat  recovery  heat  exchanger,  called  recovered 


This  balance  allows  for  the  evaluation  of  the  temperature  of  the 
outlet  gases  Tg  using  equation  (24).  The  mean  specific  heat  of  the 
outlet  gas  mixture  at  constant  pressure  is  defined  as: 


cp,s(T6 ) 


ns 


(27) 


where  j  indicates  each  chemical  compound  contained  in  the  outlet 
gas  mixture. 


Table  1 

Pressure  drop  in  system  components,  according  to  [23], 


System  component 

Pressure  drop  [kPa] 

Air  filter 

1 

Heat  exchangers  (hot  and  cold  sides) 

10 

After-burner 

2 

SOFC  (fuel  side) 

2 

SOFC  (oxidant  side) 

3 

Pipes  and  elbows 

0.5 
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thermal  power  and  represented  by  Qrec.  This  is  the  thermal  power 
that  the  system  supplies  to  an  external  user. 

The  heat  exchangers'  efficiency  can  be  expressed  as  the  ratio 
between  the  thermal  power  actually  transferred  from  hot  to  cold 
fluid  and  the  maximum  thermal  power  that  could  be  transferred: 


Cex  : 


cp.h  (7/i,i)  '  Ph.i  cp,h  ( Th.o ) '  Pfi.o 
cp,h(Th,i )  —  cp,h(Pc,i)  Tcj 


(32) 


3.5.  Power  electronics  model 

The  power  electronics'  efficiency  tjpe  is  considered  constant 
even  if  the  load  varies,  since  the  input  voltage  and  current  of  the 
equipment  are  constant,  and  the  efficiency  variation  of  DC/DC  and 
DC/AC  devices  between  20%  and  100%  of  maximum  power  is 
typically  lower  than  3%.  As  a  consequence,  it  is  possible  to  calculate 
the  net  electric  power  available  after  the  conversion  by  subtracting 
the  power  consumed  by  the  auxiliary  devices  from  the  electric 
power  generated  by  the  SOFC,  and  multiplying  the  result  by  the 
power  electronics’  efficiency: 

Pel  =  VPE'  (PeLSOFC  —  Pb  —  Pp, ethanol  ~  Pp, water)  (33) 

Using  the  net  electric  power  it  is  also  possible  to  determine  the 
global  electric  efficiency  of  the  system,  defined  as: 


Vel 


Pel 

Ptot 


(34) 


4.  Results  and  discussion 


A  computational  simulation  is  made  for  the  ensemble  SOFC-BoP 
using  the  equations  listed  in  Section  3  in  order  to  evaluate  the 
performance  of  different  potential  topologies  for  system  con¬ 
struction,  and  to  analyze  the  impact  of  the  operational  parameters 
on  system  behavior.  Since  the  number  of  parameters  is  quite  big, 
some  of  them  are  fixed  at  a  constant  value  in  order  to  limit  the 
analysis  to  a  few  relevant  situations.  The  values  chosen  for  those 
parameters  are  listed  in  Table  2. 

Besides  these  parameters,  other  variables  are  estimated  either 
by  using  data  from  the  literature,  or  by  fixing  them  at  conventional 
values.  Table  3  shows  these  variables  and  their  values. 

In  this  section,  an  indicator  called  sensibility  has  been  used  to 
investigate  the  effect  of  a  change  in  the  operational  parameters  of 
the  system  on  the  electric  efficiencies.  The  sensibility  of  the  electric 
efficiency  r]e\  to  a  change  in  the  parameter  p  from  a  value  to  a 
value  p ^  is  defined  as  the  ratio  between  the  relative  change  in  the 
efficiency  and  the  relative  change  in  p: 


4.1.  System  supplied  with  pure  ethanol 

The  stack  receives  a  flow  of  pure  ethanol  in  vapor  form  at  the 
anode  inlet.  In  the  anodic  compartments  the  fuel  is  decomposed 
into  simpler  substances,  which  take  part  in  the  electrochemical 
reactions.  Ethanol  decomposition  happens  through  pyrolysis  and 
steam  reforming.  It  is  assumed  that  the  presence  of  water  (product 
of  the  reactions),  together  with  the  materials  used  in  the  system's 
construction  and  the  structure  and  management  strategy  of  the 
BoP,  prevent  solid  carbon  deposition. 

The  main  parameters  that  affect  system  operation  are  the  A  ra¬ 
tios  for  ethanol  (Ac2h6o)  and  air  (Aajr). 


ac,h6o  = 


nC2H6Q 

ric2H6o,st 


(36) 


(37) 


The  numerator  is  the  actual  molar  flow  of  ethanol  or  air,  while 
the  denominator  is  the  minimum  molar  flow  required  to  generate 
the  desired  current  ( st  subscript). 

System  simulation  was  done  using  different  Ac2h6o  and  Aajr 
values,  obtaining  Figs.  2—4.  The  efficiency  values  are  calculated 
using  ethanol's  higher  heating  value.  Figs.  2  and  3  show  that  an 
increase  in  AC2h6o.  he.  of  fuel  excess,  causes  a  considerable  drop  in 
the  electric  efficiencies,  both  global  and  of  the  stack.  In  fact,  the 
simulations  are  operated  keeping  Pe i  constant;  therefore,  an  in¬ 
crease  in  >>c2h6o  means  that  a  higher  fuel  flow  is  being  used  to 
produce  the  same  level  of  electric  power,  resulting  in  lower  electric 
efficiencies.  The  impact  of  an  increase  in  AC2h6o  on  the  electric  ef¬ 
ficiencies  is  higher  for  low  AC2h6o  values.  Moreover,  the  effect  is 
slightly  stronger  on  the  electric  efficiency  of  the  stack  than  on  the 
global  one.  This  last  result  can  also  be  confirmed  by  calculating  the 
sensitivity  of  the  electric  efficiencies  to  the  changes  in  Ac2h6o-  In 
fact,  the  mean  value  of  the  sensitivity  of  pei  to  the  simulated 
changes  in  AC2h6o  Is  -0.426,  while  the  same  quantity  for  r)ei,S0FC  is 
equal  to  -0.448.  Both  values  are  negative,  meaning  that  an  increase 
in  Ac2h6o  leads  to  a  decrease  in  the  electric  efficiencies,  but  pei,soFC 
appears  more  sensitive  to  the  parameter  change.  Another  conse¬ 
quence  of  an  increase  of  Ac2h6o  with  Pe\  constant  is  that  a  higher 
share  of  the  fuel  energy  content  is  converted  to  heat;  therefore,  the 
thermal  efficiency  of  the  system  increases,  as  can  be  seen  in  Fig.  4. 
Also,  when  Aajr  rises,  it  is  possible  to  observe  a  limited  increase  in 
the  electric  efficiency  of  the  stack  and  a  drop  in  the  global  electric 
efficiency.  This  occurs  because  when  Aajr  rises  two  conflicting  ef¬ 
fects  develop:  the  stack’s  electric  efficiency  increases  due  to  the 
higher  partial  pressure  of  the  air  in  the  cathode;  the  blower's  power 
consumption  also  grows,  since  higher  Aa;r  implies  a  higher  air  flow. 
The  first  effect  contributes  to  an  increase  in  the  power  produced  by 


C  _  ^el 

p  pP)-p(>) 

(35) 

Table  3 

pd) 

Estimated  variable  values. 

Table  2 

Variable 

Value 

Unit 

Fixed  parameter  values. 

Power  electronics'  efficiency 

94 

m 

Operational  parameter 

Value 

Unit 

Blower's  isentropic  efficiency 

Pumps'  isentropic  efficiency 

70 

70 

m 

m 

Total  electric  power  of  the  system 

500 

[W] 

Pressure  drop  in  the  cathodic  line 

47.5 

[kPa] 

Single  cell  area 

50 

[cm2] 

Pressure  drop  in  the  anodic  line 

55 

[kPa] 

Number  of  single  cells 

30 

Ambient  pressure 

101.325 

[kPa] 

SOFC  inlet  gas  temperature 

600 

ra 

Ambient  temperature 

25 

ra 

SOFC  internal  temperature  (solid  components) 

850 

ra 

Dispersed  thermal  power  (SOFC  and  after-burner) 

120 

[W] 

System  outlet  gas  temperature 

80 

ra 

Dispersed  thermal  power  (air  pre-heater) 

80 

[W] 

Inlet  water  temperature 

25 

ra 

Dispersed  thermal  power  (fuel  pre-heater) 

80 

[W] 

Outlet  water  temperature 

45 

ra 

Dispersed  thermal  power  (heat  recovery  heat  exchanger) 

50 

[W] 
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Fig.  2.  Global  electric  efficiency  as  a  function  of  /.c,Mii0.  The  graph  shows  three 
different  trends  in  the  global  electric  efficiency,  each  one  corresponding  to  a  different 
value  of  Aajr. 


electric  efficiency.  Compared  to  AC2h6o.  ^air  influences  the  system 
performance  in  a  much  smaller  way,  especially  in  simulations 
operated  with  high  XC2h6o  values.  The  mean  value  of  the  sensitivity 
of  77ei  to  the  simulated  changes  in  Aair  is  -0.018,  meaning  that  an 
increase  in  Aajr  leads  to  a  much  smaller  decrease  in  7)e i  than  the  one 
observed  in  the  case  of  an  increase  in  AC2h6o- 

According  to  the  analysis  conducted  above,  in  order  to  maximize 
the  electric  efficiency,  it  is  necessary  to  reduce  AC2h6o  as  much  as 
possible.  The  lowest  value  of  XC2h6o  that  it  will  be  possible  to  use  in 
the  real  system  mainly  depends  on  the  actual  constructive  char¬ 
acteristics  of  the  stack.  Ac2h6o  values  lower  than  1.2  are  not 
considered  in  the  simulations,  because  below  this  threshold  the 
risk  of  a  fuel  shortage  for  the  electrochemical  reactions  would  be 
substantial.  In  fact,  the  reactions  that  happen  in  the  cells  are  never 
ideal,  and  therefore  it  is  always  necessary  to  ensure  some  fuel 
excess  to  be  sure  that  the  stack  is  able  to  sustain  the  desired  current. 
Furthermore,  a  minimal  difference  of  0.2  between  Ac2h6o  and  Aajr  is 
maintained,  in  order  to  guarantee  complete  fuel  combustion  inside 
the  after-burner.  In  fact,  a  smaller  air/fuel  ratio  could  result  in  the 
emission  of  uncombusted  hydrocarbons  from  the  system  outlet. 


4.2.  System  supplied  with  ethanol  and  water 

Solid  carbon  deposition  inside  the  SOFC  and  the  pipes  is  an 
important  problem  for  SOFC-based  electric  generation  systems.  To 
avoid  it,  the  fuel  is  often  supplied  to  the  system  with  a  certain 
amount  of  water.  The  necessary  amount  of  water  changes  with  the 
system  characteristics.  The  parameter  used  to  define  this  amount  is 
the  ratio  between  steam  and  ethanol  molar  flows,  S/E 


the  stack  for  a  constant  fuel  flow,  while  the  second  implies  higher 
power  consumption  by  the  auxiliary  devices.  However,  the  second 
effect  is  more  relevant  than  the  first,  resulting  in  a  drop  in  the  net 
electric  power  produced  by  the  system  and,  therefore,  in  the  global 


S  _  Hh2o 
E  dc2H6o 


(38) 


Values  normally  presented  in  the  literature  or  used  in  prototype 
operation  are  between  2  and  5.  Fig.  5  shows  the  effect  of  an  increase 
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Fig.  3.  Stack  electric  efficiency  as  a  function  of  >c2h6o-  The  graph  shows  three  different 
trends  in  the  stack  electric  efficiency,  each  one  corresponding  to  a  different  value  of 


Fig.  4.  Global  thermal  efficiency  as  a  function  of  Xc2h6o-  The  graph  shows  three 
different  trends  in  the  global  thermal  efficiency,  each  one  corresponding  to  a  different 
value  of  Aair. 
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in  the  S/E  value  on  both  the  global  and  the  stack's  electric 
efficiencies. 

The  graph  is  plotted  using  the  couple  of  values  of  1Vc2h6o  and  7ajr 
that  generates  the  highest  electric  efficiency  for  the  case  of  pure 
ethanol,  and  under  the  assumption  that  a  pump  is  used  to  send  the 
water  flow  to  the  stack.  Fig.  5  clearly  shows  that  the  dilution  of 
ethanol  with  water  results  in  a  drop  in  the  global  and  stack's 
electric  efficiencies.  The  drop  in  the  global  efficiency  is  mainly 
caused  by  the  drop  in  the  electric  efficiency  of  the  stack,  but  there  is 
also  a  very  limited  additional  deleterious  effect,  which  is  the  power 
consumed  by  the  pump.  The  more  diluted  the  ethanol,  the  greater 
the  loss  in  efficiency.  With  the  lowest  value  of  S/E  used,  the  dif¬ 
ference  between  efficiencies  in  the  cases  of  pure  and  diluted 
ethanol  is  found  to  be  equal  to  2.8%.  This  difference  increases  as  the 
value  of  S/E  grows.  This  is  reflected  also  by  the  values  of  the 
sensitivity  of  rjei  to  the  changes  in  S/E:  the  mean  of  the  sensitivity 
values  calculated  for  the  simulated  steps  of  S/E  is  -0.058,  meaning 
that  an  increase  in  S/E  causes  a  drop  in  rjei.  This  drop  is  smaller  than 
the  one  caused  by  an  increase  in  Xc2h6o  and  bigger  than  the  one 
caused  by  an  increase  in  7ajr,  as  can  be  noted  by  comparing  the 
sensitivities  calculated  for  XC2h6o.  ^air  and  S/E.  Moreover,  the 
sensitivity  increases  when  the  considered  S/E  values  increase, 
confirming  that  the  pejorative  effect  of  a  rise  in  S/E  gets  stronger 
when  S/E  rises. 


gases  to  the  anode  inlet  line.  Through  this  mechanism,  the  steam 
produced  in  the  electrochemical  reactions  is  partially  reintroduced 
into  the  stack,  increasing  the  water  concentration.  This  is  often 
done  using  an  ejector,  which  uses  the  pressure  of  the  fresh  fuel  flow 
(supplied  from  outside  the  stack)  to  generate  the  recirculation  flow. 
The  parameter  that  fixes  the  degree  of  anodic  recirculation  is  the 
fraction  of  the  water  flow  introduced  from  outside  the  stack,  i.e.  the 
ratio  between  the  water  flow  supplied  with  the  fresh  fuel  and  the 
total  water  flow  that  enters  the  anode  (the  sum  of  the  flow  coming 
from  outside  the  stack  and  the  flow  resulting  from  gas 
recirculation): 


/h2o  = 


nH20,out 

^H20 


(39) 


The  lower  this  ratio,  the  greater  the  recirculation  gas  flow. 
System  complexity  is  minimized  when  the  stack's  water  demand  is 
satisfied  exclusively  through  gas  recirculation,  because  in  this  case 
no  water  pump  is  necessary. 

It  is  important  to  highlight  that  part  of  the  fuel  flow  defined  by 
^c2h6o  is  contained  in  the  recirculation  gases:  in  fact,  these  gases 
contain  non-oxidized  fuel  (in  the  present  model,  solely  hydrogen, 
generated  through  ethanol  decomposition  in  anodic  compart¬ 
ments)  that  can  still  take  part  in  the  electrochemical  reactions.  This 


means  that  when 

using  anodic  recirculation  an  increase  in  fuel 

4.3.  Anodic  recirculation 

One  method  that  is  used  to  increase  the  water  concentration  in 

S/E=2 

the  anode  compartments,  reducing  the  probability  of  solid  carbon 
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Fig.  5.  Electric  efficiencies,  ethanol-water  mixture.  The  trends  of  the  stack  electric 
efficiency  and  of  the  global  electric  efficiency  as  a  function  of  S/E  are  reported  for  the 
case  in  which  the  system  is  fueled  with  an  ethanol-water  mixture,  without  the  use  of 
anodic  recirculation. 


Fig.  6.  Electric  efficiencies,  anodic  recirculation,  S/E  =  2.  The  trends  of  stack's  electric 
efficiency  and  of  global  electric  efficiency  as  a  function  of/n2o  are  reported  for  the  case 
in  which  the  system  is  fueled  with  an  ethanol-water  mixture,  using  anodic  recircu¬ 
lation  and  for  a  S/E  value  equal  to  2. 
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excess  is  related  to  an  increase  in  nC2H6o  that  is  smaller  than  the  one 
described  in  Subsections  4.1  and  4.2. 

Figs.  6  and  7  show  the  electric  efficiencies  of  the  stack  and  the 
SOFC,  for  different  values  of  /H2q  and  for  S/£  equal  to  2  and  4. 
Comparing  this  data  with  the  results  displayed  in  Subsections  4.1 
and  4.2,  it  can  be  concluded  that  with  the  use  of  anodic  recircula¬ 
tion  it  is  possible  to  obtain  a  high  global  electric  efficiency,  while 
maintaining  S/£  at  the  desired  value.  When  S/£  is  equal  to  2  (Fig.  6) 
and  /h2o  is  lower  than  0.2,  the  global  electric  efficiency  is  higher 
than  with  pure  ethanol  (Fig.  2),  even  though  the  difference  is  quite 
small  (0.4%  at  the  maximum). 

As  demonstrated  in  Subsection  4.2,  an  increase  in  S/£  causes  a 
global  electric  efficiency  drop:  when  S/£  is  equal  to  4  (Fig.  7),  the 
global  electric  efficiency  is  the  same  as  that  observed  for  pure 
ethanol  (Fig.  2)  only  when/H2o  is  0.  This  means  that  for  higher  S/£ 
values  it  is  impossible  to  avoid  a  global  electric  efficiency  drop  with 
respect  to  the  pure  ethanol  case,  even  using  anodic  recirculation. 

By  calculating  the  sensitivities  of  jjei  to  the  steps  of/H2o  between 
the  various  simulations,  it  is  possible  to  retrieve  some  additional 
information.  The  average  of  the  sensitivities  is  -0.039  in  the  case  of 
S/£  equal  to  2,  while  the  same  value  is  -0.054  for  S/£  equal  to  4.  This 
suggests  that  the  influence  of  /H2q  on  r/e\  is  bigger  for  high  S/£ 
values.  Moreover,  the  absolute  value  of  the  sensitivity  is  higher 
when  the  considered  values  of/H2o  are  high.  Therefore,  the  positive 
effects  of  anodic  recirculation  are  very  evident  even  with  only  small 
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reductions  in  the  external  water  flow  (i.e./H2o  close  to  1),  while  the 
additional  benefits  are  reduced  when  /h2o  approaches  0  (at  which 
point  water  is  completely  supplied  through  anodic  recirculation).  In 
this  case,  as  in  Subsection  4.2,  the  sensitivity  of  7)ei  to  changes  in 
/h2o  proves  lower  than  the  sensitivity  of  jjei  to  changes  in  XC2h6o- 


4.4.  Cathodic  recirculation 

As  with  the  anode,  it  is  possible  to  recirculate  part  of  the  cathode 
outlet  gases  to  the  inlet  of  the  cathodic  compartments.  In  this  case 
an  ejector  is  often  used.  The  cathodic  recirculation  reduces  the 
compressed  air  flow  required  to  obtain  a  given  value  of  Aajr:  this 
reduces  the  thermal  energy  required  to  pre-heat  the  air  and  the 
power  consumption  of  the  blower.  As  a  consequence,  when  using 
this  solution  the  auxiliary  devices  power  consumption  drops  and 
the  system’s  thermal  performance  improves. 

In  this  case,  the  parameter  that  characterizes  system  operation 
is  the  cathodic  recirculation  fraction,  defined  as  the  ratio  between 
recirculated  gas  flow  and  total  cathode  outlet  flow. 


Irec.c 


tirec.C 

tiout.C 


(40) 


An  increase  of  /REc,c  results  in  a  reduction  of  oxygen  concen¬ 
tration  in  the  gas  mixture  that  enters  the  cathodes,  while  the  ni¬ 
trogen  content  of  the  same  mixture  increases.  Therefore,  /rec.c  can 
assume  only  those  values  that  allow  for  a  minimum  flow  of  oxygen 
to  sustain  the  reactions  in  the  SOFC  and  in  the  after-burner. 

Fig.  8  shows  the  effect  of  changes  in  /rec,g  as  expected,  the  in¬ 
crease  in/REc,c  causes  a  decrease  in  blower  power  consumption  and 
an  increase  in  recovered  thermal  power  (due  to  a  decrease  in  the 
thermal  power  required  to  pre-heat  the  air).  The  cathodic  recircu¬ 
lation  therefore  reduces  the  electric  consumption  of  the  auxiliary 
devices,  while  improving  the  thermal  management  of  the  system 
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Fig.  7.  Electric  efficiencies,  anodic  recirculation  S/E  =  4.  The  trends  of  stack's  electric 
efficiency  and  of  global  electric  efficiency  as  a  function  of/n2o  are  reported  for  the  case 
in  which  the  system  is  fueled  with  an  ethanol-water  mixture,  using  anodic  recircu¬ 
lation  and  for  an  S/E  value  equal  to  4. 


Fig.  8.  Electric  power  consumed  by  the  blower  and  recovered  thermal  power  using 
cathodic  recirculation.  For  the  case  in  which  the  system  employs  the  cathodic  recir¬ 
culation,  trends  of  the  electric  power  consumed  by  the  blower  and  the  recovered 
thermal  power  as  a  function  of  /rec,c  are  reported. 
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Fig.  9.  Cathodic  recirculation:  stack  electric  efficiency,  global  electric  efficiency  and  global  thermal  efficiency.  The  trends  of  stack's  electric  efficiency,  of  global  electric  efficiency  and 
of  global  thermal  efficiency  as  a  function  of  /RE c.c  are  shown  for  the  case  of  cathodic  recirculation. 


and  consequently  the  global  thermal  efficiency.  Fig.  9  shows  the 
behavior  of  the  stack's  electric  efficiency,  of  the  global  electric  ef¬ 
ficiency  and  of  the  global  thermal  efficiency  using  cathodic 
recirculation. 

An  increase  in  /rec.c  causes  a  reduction  in  cathodic  oxygen 
concentration  and  consequently  a  drop  in  stack  electric  efficiency. 
At  the  same  time,  there  is  a  reduction  in  the  energy  consumption  of 
the  auxiliary  devices.  Both  effects  influence  the  global  electric  ef¬ 
ficiency,  but  the  drop  in  the  electric  efficiency  of  the  stack  domi¬ 
nates,  resulting  in  a  global  electric  efficiency  drop  when  /rec,c 
increases.  However,  this  drop  in  global  electric  efficiency  is  very 
small,  around  0.2%  change  when  comparing  the  cases  of  /rec,c 


respectively  equal  to  0  and  0.3.  The  size  of  the  impact  of/REc,c  on  i 
is  so  small  that  the  corresponding  sensitivities  are  almost  0  in  every 
case.  On  the  other  hand,  the  effect  of  cathodic  recirculation  on  the 
global  thermal  efficiency  is  stronger:  the  thermal  performance  of 
the  systems  increases  by  approximately  0.7%  for  the  same  change  in 
/rec.c- 

These  results  show  that  cathodic  recirculation  is  not  appropriate 
for  a  system  that  prioritizes  electric  production,  because  it  causes  a 
global  electric  efficiency  reduction.  However,  cathodic  recirculation 
can  be  advantageous  if  global  thermal  efficiency  is  an  important 
parameter,  because  it  generates  an  increase  in  thermal  efficiency 
that  is  higher  than  the  decrease  in  global  electric  efficiency.  Still, 


Table  4 

Comparison  between  best  performing  operating  conditions,  identified  for  each  analyzed  system  configuration. 


Variable 

1 

2 

3 

4 

Unit 

Pure 

ethanol 

Ethanol  and 
water  mixture 

Pure  ethanol  with 
anodic  recirculation 

Pure  ethanol  with 
cathodic  recirculation 

kc2H60 

1.2 

1.2 

1.2 

1.2 

7air 

1.6 

1.6 

1.6 

1.6 

S/E 

0 

2 

2 

0 

/rec.c 

0 

0 

0 

0.3 

/h2o 

1 

1 

0 

1 

Global  electric  efficiency 

43.9 

41.1 

44.3 

43.7 

i%] 

Global  thermal  efficiency 

13 

17.1 

12 

13.7 

i%] 

Ethanol  flow  rate 

2.9 

3.1 

2.87 

2.92 

[ml  min  *] 

Fuel  chemical  power 

1140 

1215 

1128 

1144 

[W] 

SOFC  DC  electric  power 

555 

557 

558 

553 

[W] 

Recovered  thermal  power 

149 

208 

135 

156 

[W] 

Power  consumed  by  the  ethanol  pump 

4 

4 

4 

4 

[mW] 

Power  consumed  by  the  water  pump 

0 

3 

0 

0 

[mW] 

Power  consumed  by  the  blower 

23.2 

24.7 

25.7 

20.6 

[W] 

Mean  cell  voltage 

694 

653 

651 

688 

[mV] 

Stack  voltage 

20.8 

19.6 

19.5 

20.7 

[V] 

Stack  current 

26.7 

28.4 

28.5 

26.8 

[A] 

After-burner  outlet  gas  temperature 

1026 

1027 

984 

1067 

rc] 

The  global  electric  efficiencies  reported  in  bold  characters  are  related  to  the  system  configurations  indicated  as  most  convenient  in  Section  5. 
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Case 

Global  electric  efficiency 

Global  thermal  efficiency 

Pure  ethanol 

0.439 

0.13 

Ethanol+H20 

0.411 

0.171 

An  rec 

0.443 

0.12 

Cat  rec 

0.437 

0.137 

0  0.1  0.2  0.3  0.4  0.5 

Efficiency 

Global  electric  efficiency  ■  Global  thermal  efficiency 


Fig.  10.  Global  system  efficiencies,  comparison.  The  highest  values  of  global  electric 
efficiency  and  global  thermal  efficiency  obtained  for  each  system  configuration  are 
reported. 

this  increase  in  global  thermal  efficiency  is  not  very  high;  therefore, 
this  beneficial  effect  would  hardly  legitimate  the  higher  system 
complexity  required  for  the  application  of  cathodic  recirculation. 

4.5.  Comparison 

Table  4  shows  the  values  of  the  variables  analyzed,  taking  into 
account  the  different  operational  conditions  considered  in  this 
study.  The  situations  included  in  Table  4  and  Figs.  10  and  11  are  the 
ones  that  present  the  highest  global  electric  efficiency  for  each 
configuration,  except  for  the  case  of  cathodic  recirculation,  where 
the  situation  of  maximum  total  efficiency  (sum  of  thermal  and 
electric  efficiencies)  was  reported.  In  Fig.  10  it  can  be  seen  that  the 
global  electric  efficiency  shows  its  maximum  value  when  anodic 
recirculation  is  adopted,  and  that  the  variation  of  this  parameter 


among  all  situations  is  quite  small.  The  figure  also  shows  an  inverse 
relationship  between  electric  efficiency  and  thermal  efficiency.  The 
cases  characterized  by  higher  global  thermal  efficiency  also  present 
higher  values  of  recovered  thermal  power,  as  shown  in  Fig.  11.  This 
figure  also  reports  data  about  blower  power  consumption,  which  is 
approximately  constant  in  all  situations,  except  in  case  of  cathodic 
recirculation,  when  the  value  is  lower.  Also  stack's  electric  power 
varies  in  a  limited  manner  among  the  reported  cases. 

The  analysis  of  the  system  behavior  when  the  various  parame¬ 
ters  change  provides  information  about  their  desirable  values.  In 
general,  in  order  to  maximize  the  global  electric  efficiency  of  the 
system  it  is  good  to  work  with  low  >>c2h6o  and  -tar  values  and  to  set 
/rec.c  equal  to  0  (i.e.  not  to  use  cathodic  recirculation).  If  the  char¬ 
acteristics  of  the  system  impose  a  water  flow  to  the  anode  together 
with  the  fuel  to  avoid  solid  carbon  deposition  (i.e.  a  S/E  value  higher 
than  0),  it  is  desirable  to  use  values  of  S/E  and  fHl0  as  low  as 
possible.  Comparing  the  average  sensitivities  of  ?je l  to  a  change  in 
the  various  parameters,  it  is  possible  to  conclude  that  the  most 
influential  parameter  is  >-c2h6o.  followed  by  S/E  and  /h2o.  this  last 
one  especially  in  the  case  of  high  S/E  values.  Aajr  and  /rec.c  are 
characterized  instead  by  a  limited  effect,  almost  inexistent  in  the 
case  of/REC,c. 

5.  Conclusions 

The  analysis  of  the  simulation  results  with  pure  ethanol,  shown 
in  Figs.  2—4  and  reported  in  Table  4,  leads  to  the  conclusion  that  an 
increase  in  fuel  excess,  i.e.  a  higher  value  of  XC2h6o.  causes  a  sub¬ 
stantial  reduction  in  the  global  electric  efficiency  and  an  increase  in 
the  thermal  efficiency.  Fig.  5  clearly  shows  that  ethanol  dilution 
with  water  causes  a  considerable  drop  in  system  efficiency  if  water 
is  supplied  by  a  dedicated  pump.  The  effect  on  the  efficiency  of  the 
system  is  less  negative,  or  even  positive,  if  water  concentration  is 
increased  by  partial  recirculation  of  anode  outlet  gases.  The  recir¬ 
culation  of  cathode  outlet  gases  has  a  more  limited  effect,  which  is 
therefore  not  interesting  for  the  considered  application. 

Table  4  summarizes  the  values  of  the  variables  that  were  ob¬ 
tained  during  the  simulations.  All  the  efficiency  values  reported  in 
the  table  are  related  to  ethanol's  higher  heating  value.  The  results 
show  that  there  are  two  convenient  system  configurations: 

•  Operation  with  pure  ethanol,  without  recirculation 


Case 

Stack  electric  power 

Power  consumed  by  the  blower 

Recovered  thermal  power 

Pure  ethanol 

555.072 

23.154 

148.77 

Ethanol+H20 

556.603 

24.682 

208.031 

An  rec 

557.637 

25.718 

134.919 

Cat  rec 

552.545 

20.625 

156.233 

Cathodic  recircula 

tion 

Anodic  recirculati 

)n 

■ 

Ethanol-water  m 

ixture 

■ 

|  Pure  ethanol 

0 

200  400 

Power  [W| 

Stack  electric  power 

■  Power  consumed  by  the  blower 

■  Recovered  thermal  power 

Fig.  11.  System  power,  comparison.  The  values  of  stack's  electric  power,  blower  power  consumption  and  recovered  thermal  power,  obtained  for  the  same  system  configurations 
considered  in  Fig.  10,  are  reported. 
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In  this  case  the  system  is  very  simple,  since  it  does  not  require 
ejectors  or  a  water  pump.  At  the  same  time,  the  global  electric  ef¬ 
ficiency  is  high  (43.9%).  However,  system  operation  in  this  case 
requires  the  use  of  suitable  materials  and  an  appropriate  BoP 
structure  and  management,  to  prevent  solid  carbon  deposition. 

•  Operation  with  partial  recirculation  of  anodic  reaction  products 

In  this  configuration  the  system  is  more  complex,  even  if  it  does 
not  require  a  water  pump,  but  there  is  a  lower  probability  of  solid 
carbon  deposition.  Moreover,  the  global  electric  efficiency  is 
slightly  higher  than  the  one  observed  when  operating  with  pure 
ethanol  without  recirculation  (44.3%). 

Finally,  these  results  suggest  that  the  most  convenient  solution 
for  this  project  is  a  pure  ethanol  operation  without  recirculation, 
made  possible  by  the  use  of  proprietary  anode  components  ([15,16], 
),  which  were  developed  especially  to  avoid  solid  carbon  formation. 
As  a  consequence,  the  data  from  this  pure  ethanol  simulation  will 
be  used  as  a  reference  to  fix  a  tentative  operating  point  for  the 
system  and  to  continue  with  the  prototype  project,  as  established  in 
the  research  plan.  In  the  continuation  of  the  prototype  develop¬ 
ment,  economic  considerations  about  the  cost  of  the  system  will  be 
included  in  the  decision  process. 

As  stated  in  the  introduction,  the  efficiency  of  fuel  cell  based 
electricity  generation  systems  is  approximately  constant  for  a  wide 
range  of  system's  nominal  power.  This  means  that,  in  the  future,  it 
would  be  possible  to  easily  extend  these  results  to  other  power 
levels,  implementing  electric  generation  systems  characterized  by 
higher  powers. 
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